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Abstract 

The J/^(1S) and ^'(2S) production near threshold in antiproton-nucleus reactions is calculated 
on the basis of the Glauber model. The model takes into account the antiproton (pre) absorption, 
proton Fermi motion, and charmonium formation length. We confirm earlier prediction that the 
charmonium production in pA collisions at pi a b = 3—10 GeV/c is not influenced by formation 
length effects and is very well suited to determine the genuine charmonium-nucleon dissociation 
cross sections. The comparison is performed with the J/^ photoproduction at high energies, where 
formation length effects play a very important role. However, we demonstrate that the detailed 
structure of the proton and neutron density profiles has to be taken into account, if one wants to 
extract information on the J/^N dissociation cross section from J/*f> transparency ratios. These 
studies are relevant for the upcoming PANDA experiment at FAIR. 
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1. INTRODUCTION 



The exploration of the properties of quantum chromodynamics (QCD) matter under 
extreme conditions is one of the most challenging endeavors of todays high-energy physics. 
The dynamics of heavy quarks, i.e. c- and 6-quarks, and their bound states has been shown 
to provide new insights into these questions. A better understanding of the interactions of 
charmonia with cold nuclear matter is especially very important for the studies of charm 
production in heavy-ion collisions, in particular, of the J/^> suppression in a quark-gluon 
plasma The interpretation of J/ty yields from p-A and non-central AA collisions at 
moderate energies (\/s =17.3 GeV) in terms of a usual hadronic scenario requires the J/^> 



dissociation cross section on a nucleon to be about 6 — 7 mb {2, 3|. Moreover, as stated in 
Ref. [2I, such cross sections are consistent with the world data on J/ty transparency ratios 
from other reactions induced by elementary particles (7, ir, p) on nuclei. However, if one 
takes into account contribution of J/^B produced in the decays of higher charmonium states 
and larger cross section of inelastic \cN interactions, the Jlj&N dissociation cross section 
turns out to be about two times smaller, i.e. about 3.5 mb 4]. 

A problematic feature of all existing experimental data on J/ty transparency ratios, 

bA = x; ' t ) 

(here "p" denotes any kind of elementary projectile and N stands for the nucleon) in elemen- 
tary particle-induced reactions is that the charmonia are produced at high momenta from 
about 20 GeV/c up to several TeV/c. Due to the large Lorentz-boost of the charmonia, the 
existing data only give access to the interactions of prehadronic cc and ccg configurations 
with the nuclear medium which would transform into J/^> well after the target. Thus, the 
extracted value of the dissociation cross section needs to be studied in other kinematical 
situations, when the charmonium moves more slowly through the nuclear target. 

Brodsky and Mueller 5| proposed to measure color transparency effects for J/^> produc- 
tion in p-nucleus interactions, where the threshold beam momentum (pthx = 4.07 GeV/c for 
J/ty production in pp interactions) is quite low. The first theoretical study of this reaction 



has been done in Ref. 61 within the color diffusion and Glauber models tal 

r, 

the Fermi motion of the nucleons. The most important results of Ref. [6| are (i) the high 
sensitivity of the production cross section (JpA-+R{A-i)* of a charmonium state R, where R 



king into account 



stands for J/^, or Xc on a nucleus to the .RiV-dissociation cross section; 
(ii) due to the Fermi motion the cross section of charmonium production on a nucleus is 
strongly suppressed, i.e. (TpA-*R(A-\)* ~ (10~ 4 — 10~ 3 )Zapp^.R at the beam momentum of 
the on-shell R production..; 

and (hi) crpA-+R(A-i)* is sensitive to the color transparency effect for the incoming antiproton 
(which can, however, be reformulated as a phenomenological treatment of p absorption), but 
insensitive to the color transparency for the produced J/W and having modest sensitivity 
for heavier charmonium state R. 

Later-on, in work the J/^f and A c production in pA collisions at the ^' and J/^f 
production thresholds has been addressed with a focus on the effects of the non-diagonal 
transitions WN — > J/^N. However, Fermi motion effects have been neglected in Ref. 
We would also like to mention a very inspiring feasibility study for Fermi National Accel- 
erator Laboratory (Fermilab) 8j, which for the first time presented the beam momentum 
dependence of the charmonium production cross sections in pA collisions with heavy gas 
targets (CH 4 , N, O, Ne, Ar and Xe). 

The main purpose of the present work is to perform a detailed theoretical analysis of the 
J /ty production in pA collisions near threshold. This is one of the subjects of the planned 
PANDA experiment at FAIR 9|. Comparisons of the J/^> transparency ratios in pA and 
'-/A reactions are performed with a focus on the sensitivity to the J/^/N dissociation cross 
section. 

Section [2] explains our model. The model predictions for the charmonium production in 
pA reactions are given in Sec. [3] with an emphasis on the sensitivity to the charmonium- 
nucleon dissociation cross sections. For comparison, we have also calculated the J/ty trans- 
parency ratios in photo-induced reactions and showed that the existing experimental data 
from Stanford Linear Accelerator Center (SLAC) for E-y = 20 GeV 10j and, moreover, of 



Fermilab for £ 7 = 120 GeV 



11] do not allow to constrain the J/^fN dissociation cross 



section due to the formation length effects and large experimental errors. Finally, in order 
to evaluate the influence of multistep processes on J/^> production in pA reactions the re- 
sults obtained within the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) model [3] are 
presented. The conclusions are given in Sec. HI 

1 This estimate corresponds to the perfect beam resolution for pp collisions. See Eq. ([20|) below and the text 
after it. 
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2. MODEL 



In the reaction pA — > R + X at beam momentum close to the charmonium R produc- 
tion threshold, the produced charmonium carries nearly the entire antiproton momentum. 
Therefore, we can apply a Glauber model similar to that of Ref. js]. The cross section of 
the charmonium R production in a pA collision is given by 

oo oo 

0pA-+R(A-i)* = 2vr J dbbvp 1 J dzPp iSUVV (z,b)T p ^ R (z,b)P RiSurv (z,b) , (2) 

-oo 

where the integration is done over the antiproton impact parameter b and the longitudinal 
coordinate z. v p = p\ab/E p is the antiproton velocity with respect to the target nucleus. The 
in-medium width of the antiproton with respect to the charmonium production is 

f 2d p 

Tp^ R (z,b) = J j^-^Vp p ap p ^ R x(p,Pp)fp(z,b,p) , (3) 



where Vpp = q^fs j E p E p is the antiproton-proton relative velocity with q = y s/4 — m? be- 
ing the center-of-mass (cm.) momentum of the antiproton and the proton at the cm. 
energy of a/s; m = 0.938 GeV is the nucleon mass. To simplify the discussion we use 
here the simplest model for the proton momentum distribution, i.e. the Fermi distribution 
f p (z,b,p) = Q(pf, p — |p|)) where pp, p = (3n 2 p p (z, 6)) 1//3 is the local Fermi momentum of 
protons, and p p is the local proton density. The more realistic spectral function with corre- 
lations would slightly reduce the value of cross section near maximum and add momentum 
tails removing the sharp cutoffs. Currently the program of the experimental studies of the 
short-range correlations is under way at the TJNAF and models of the nuclear spectral 



functions incorporating these findings are being developed, for a recent review see [13| . This 
would allow in the near future to perform more accurate calculations of the rates of the 
charm production for the kinematics where p produces J/^f in the interaction with a fast 
nucleon. 

If the beam momentum is close to that of exclusive R production at the mass pole, then 
all processes except pp — > R can be neglected, and one can use in Eq.Q the exclusive 
resonance production cross section <J PP ^. R instead of the inclusive one, a pp ^ R x ■ For a pp ^ R 
we apply the relativistic Breit-Wigner formula, 

3tt 2 



a 



Pp^R 2 



■y/sT R ^ pp A R (s) , (4) 



with the resonance spectral function 



1 yfiT R 
it (a- mlf + a I* 



The remaining two important ingredients of Eq.(j2]) are the survival probability of the 
antiproton until it reaches the point (z,b), 

Pp,smv{z,b) =expj- J dz'p(z',b)a^\pi ah )^ , (6) 

and the survival probability of the charmonium R until it is emitted to the vacuum 

PR,swrv(z, b) = exp j- J dz'p(z', b)a RN (p R , z' - z) j . (7) 

In Eqs.fjSJ),©, p = p P + p n is the total nucleon density, cr^ 1 is the pN inelastic cross section 

a pN (Plab) = ~ a pN 5 ( 8 ) 

and a R s N is the charmonium-nucleon effective cross section, which will be explained be- 
low. The total and elastic antiproton- neutron cross sections in Eq.(jHJ) are set equal to the 
antiproton-proton ones. The pp cross sections are taken from the PDG parametrization jl^ : 

*j?(Pbb) = 38.4 + 77.6pf a °- 64 + 0.26 ln 2 (p lab ) - 1.2 ln(p lab ) , (9) 
<(piab) = 10.2 + 52.7pf ab - 16 + 0.125 In 2 (p lab ) - 1.281n(p lab ) , (10) 

where the beam momentum, pi ab , is in GeV/c and the cross sections are in mb. 

Let us next turn to the time dependence of the charmonium formation. This is expressed 
via the charmonium-nucleon effective cross section, a R s N (p R , z), which is a function of the 
charmonium momentum p R in the target nucleus rest-frame and of the distance z from the 



cc-pair production point. Following Refs. [6j, ll5| we express a e RN in terms of a formation 



length l R : 



z \ T < n 2 kf > ( / z 



l R J m% V VZr, 

with r — 1. In Eq.flTTJ, n is the number of hard gluons in the intermediate state and 
< k 2 t > 1 / 2 - 0.35 GeV/c is the average transverse momentum of a quark in a hadron. 
Assuming that the reaction pp — > R is dominated by qqq + qqq annihilation into three 



®{Ir-z) + Q{z-Ir)\ (11) 



hard gluons jsj, we will use the value n — 3. The formation lengths of hadrons are model- 
dependent. For the J/\l/ we apply a standard formula with an energy denominator (c.f. 



113) 

For we rely on the estimate of Ref. {4] 



- ^ 2 2pj/ : 2 • ( i2 ) 



which, however, has large theoretical uncertainty. 

Formula ([3]) for the partial width of the antiproton with respect to the process pp —> R 
can be simplified in the limit of small width of the resonance R. To this aim we perform the 
integration over proton momentum in ([3]) using the spherical coordinate system with z-axis 
along the antiproton beam momentum: 

PF,p 1 

IV* = / dpp 2 J d cose ^A R (s) , (14) 

0-1 ^ 

where cos = p z /p, and s = (E p + .Ep) 2 — p 2 — p 2 ah — 2pp lab cos 0. It is convenient to make 
the angular dependence in the spectral function explicit: 

1 

Ar{s) = 47r mab [(A(p)-cos0) 2 + 7 2 /4] ' (15) 
where A(p) = [(E p + E p ) 2 - p 2 - p 2 Xab - m 2 R ]/2pp iah and 7 = y/sT R /pp lab - If 7 < 1, one can 
replace the Breit-Wigner distribution by the 5-functional distribution □ 

A R (s) ~ -J— «J(A(p) - cos 0) , (16) 

and set ~ m R in Eq. (fT4j) . Substituting Eq. ljTB"]) in Eq. (I14[) and performing the integra- 
tion over cos we obtain the following formula: 

min(p2,PF,p) 



mm 



in(pi,p F ,p) 

where q R = \Jm 2 R /A — m 2 . The limiting momenta, p\ and P2, are, respectively, the smaller 
and larger solutions of the equation A(p) = ±1. 



2 Since p ~ pi? p ~ 0.3 GcV/c, we obtain the estimate 7 ~ 10 -4 for the J/\f and charmonium states. 
We have also checked in some selected cases that the direct (however, extremely CPU-time consuming) 
Monte-Carlo calculation of the momentum integral in Eq.Q gives results indistinguishable from those 
obtained assuming zero width of charmonium states. 
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In order to proceed further, we have to specify the dispersion relation between the energy 
and momentum of a proton in the target nucleus. The choice consistent with a model 
where nucleons carry all nucleus momentum in the infinite momentum frame is to set the 
proton energy constant independent on the proton momentum, i.e. E p = m — B , where 
B ~ 8 MeV is the nucleus binding energy per nucleon (in actual calculations we used the 



nucleus-dependent empirical values of the binding energies from Ref. [16|). In this case, 
using expression v pp = q R m R / E P E P allows us to take the momentum integral in Eq.(Tj7 
analytically: 



P ^ R %P\abEpE p q R 



mm(p 2 ,p F:P ) 2 - min(pi,p FiP ) 2 ] , (18) 

with pi t 2 = \piahT \J (E p + E p ) 2 — m 2 R \. Since p 2 3> Pf, p , we can replace the upper integration 
limit in Eqs . f 1 1 7 j) . ( j 1 8 j) by Pf,p- At the beam momentum of the on-shell R production on the 
proton in vacuum at rest, p\ ~ and Eq. ffTgj) simplifies to 

pon-sheU ^ 3m R T R ^ Pp p Fp 
P ^ R %V\abEpE p q R 

Thus r^^ he11 oc p 2 / 3 . The deviation from the usual linear density dependence originates 
from the narrowness of the resonance state R: Owing to the Fermi motion it is difficult to 
find a proton which exactly matches the on-shell resonance kinematics. 



Equation ( fl9l) leads to the estimate of Ref. js] for the ratio 

(T P A-+R(A-i)* _ 3irm R mT R ^ w _ 4 
Zcrpp^R ~ 4(m| - 2m 2 )v p p F , p 

for r# ~ 93 keV in the case of J/^>. Such a strong reduction implies, however, the antiproton 

energy being precisely on the R on-shell peak, i.e. E p = m 2 R /2m — m (or p^ = m R q R /m). 

If the beam energy resolution AE does not allow to resolve the on-shell R production 

in the pp — )• R reaction, i.e. AE ^> T R , the r.h.s. of Eq. ([20]) should be multiplied by 

2mAE/irm R T R (see also Ref. {fl). We emphasize that the beam energy resolution is strongly 

influencing the elementary cross section <J PP -> R , but not the cross section on the nucleus 

a p A->R(A-i)* , since the latter changes on the rather large scale AE ~ pf, p ~ 0.3 GeV only. 

3. RESULTS 

We will consider the following target nuclei: 9 Be, 12 C, 16 0, 27 A1, 40 Ca, 56 Fe, 63 Cu, 75 As, 
ii2,ii6,i20,i24 Sll) i42 Ce; isi^ 197 Au and 208p b> Thig choice ig mostly mot i va t e d by the avail- 
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ability of neutron density parameters 17H19I]. 

For light nuclei (A < 20) we use the proton and neutron density profiles of the harmonic 
oscillator model 

2" 

exp{-(r/R q ) 2 } , q 



1 + a 



/ r 

q I ~~R 

tin 



p, n 



(21) 



For heavy nuclei (A > 20) we apply the two-parameter Fermi distributions 

-l 



Pq{r) = p\ 



exp 



R„ 



1 



q = p,n 



The normalization constants p° q are chosen such that 

d\ Pp {r) = Z , J d 3 rp n {r) — A — Z . 
The charge-density-distribution parameters are taken from a standard compilation 



(22) 



(23) 



20|. The 



neutron density parameters for most of nuclei are taken from Nieves et al. [18j , who report 
the fits to the Hartree-Fock calculations with the density-matrix expansion 21] . In the 
present calculations we use the point density parameters of protons and neutrons which 
were obtained from charge density parameters and neutron matter density parameters by 
employing the correction formulas from Ref. 18] . 

The neutron density parameters for some of nuclei used in our calculations are, however, 
not given in Ref. [18|. For 9 Be and 181 Ta we rely upon the Glauber model analysis of 1 GeV 
proton charge-exchange scattering by Koptev et al fl7|. while for Sn isotopes we employ 
the results of the antiprotonic X-ray analysis of Schmidt et al 19[. The neutron density 
parameters of these nuclei are collected in Table H] 

The main features of the antiproton-nucleus interaction with a heavy nucleus leading to 
the exclusive J/^> production are visualized in Fig. [lj which shows the p survival probability 
-Pp,surv; the partial J/^> production width T^j/^ and their product as functions of z at the 
two different values of an impact parameter for the 181 Ta target. The beam momentum 
4.07 GeV/c is chosen to set the produced J/^> on-shell for the proton target at rest. Thus, 
according to Eq. (lT91) . Tp^j/^ oc p 2 J z . As expected, the antiproton is almost completely 
absorbed in the diffuse surface region, where the partial width Tp^j/^, is relatively small. 
Therefore, p-absorption strongly (~ 5 times) reduces the J/^f production governed by the 
product Pp,sur v rp-!>j/*- Moreover, the surface absorption of p leads to the significant sensi- 
tivity of J/ty production to the diffuseness of the proton density distribution, in-particular, 



TABLE I. Neutron density parameters (in fm) for some of the nuclei used in calculations. For 
9 Be and 181 Ta nuclei, the neutron matter density while for Sn isotopes the point neutron density 
parameters are given. 

Nucleus R n a n 

9 Be 2.11 1.000 

181 Ta 6.42 0.640 

112 Sn 5.416 0.543 

116 Sn 5.399 0.552 

120 Sn 5.356 0.565 

124 Sn 5.530 0.558 



for peripheral collisions (cf. thick and thin lines in the lower panel of Fig. [I]). This sensitivity 
is more clearly demonstrated in Fig. [2j where the impact parameter dependence of the J/^> 
production probability (neglecting J/\l/ absorption) 

oo 

Pj/*(6) = vp 1 J dzP p 

,surv (z, b)Tp^j/i S ,(z, b) (24) 

— oo 

is shown for the two slightly different values of the charge density diffuseness parameter. 
The same effect shows up also in the mass dependence of the J/^> transparency ratio (upper 
panel of Fig. [4] below). Let us now discuss the impact-parameter integrated cross sections. 

All cross sections have been calculated assuming the minimum bias triggering condition 
for the p-nucleus collisions. Numerically, this has being done by setting the upper limit for 
the impact parameter integration in Eq. (J2J) equal to a large value, R n + 10a n for all nuclei 
except 9 Be, and 8 fm for 9 Be. 

Fig. [3] shows the beam momentum dependence of J/ty production cross section for several 
target nuclei. Apart from the result without J/^> absorption, we present two calculations 
with different J/tyN dissociation cross sections. The choice <Jj/^n — 3.5 mb is motivated by 
the early experiment on J/\P photoproduction at E 1 = 20 GeV at SLAC [10], while the value 
aj/iSN ~ 6 mb is obtained in Q from the global Glauber fit of the J/\I/ transparency ratios 
in high-energy 7-, p-, p- and 7r-induced reactions. The large J/^/N inelastic cross section in 
the range 6 — 8 mb is reported in recent calculations employing effective Lagrangians of the 



local hidden gauge theory 



22|. Since the J/\l/-formation length ~ 0.4 fm at p lab = 4 
9 



p(4.07 GeV/c) 181 Ta -> J/* (A-1) 
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FIG. 1. (Color online) Dependence of the antiproton survival probability (Eq.(|6|), solid lines), of 
the J/ty production width (Eq. (|19p . dash-dotted lines) and of their product (dashed lines) on the 
longitudinal coordinate z for the central (b = 0.62 fm, upper panel) and peripheral (b = 6.8 fm, 
lower panel) collisions p+ 181 Ta at pi a b = 4.07 GeV/c. The width is given in units of 10 _8 c/fm. 
The thick lines are obtained with the charge density diffuseness parameter a c h = 0.64 fm, while 
the thin lines - with a c h = 0.52 fm. The center of the nucleus is at b = 0, z = 0. The antiproton 
propagates in the positive z direction. 

GeV/c, the results are practically insensitive to the formation length effects, and we show 
only calculations with Ij/y = in Eq.f lllj) . On the other hand, the final J/ty yield reveals 
a clear sensitivity to the J/^/N dissociation cross section which becomes more pronounced 
for heavier targets. 
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FIG. 2. (Color online) J/^f production probability (Eq. (l24p ) multiplied by impact parameter b as 
a function of b for p+ 181 Ta collisions at pi a b = 4.07 GeV/c. The lines with solid (open) circles are 
calculated with the charge density diffuseness parameter a c h = 0.64 (0.52) fm. J absorption is 
turned off. 



In Fig. S]we show the transparency ratio, 



CpA^R(A-l)* /27\ 2 / 3 



A) 



(25) 



calculated with the J/\l/ production cross sections at their peak values (see Fig. [3]) and 
rescaled by A~ 2 ^ 3 . This rescaling factor corresponds to the surface-dominated p absorption 
at moderate beam momenta. The nucleus 27 A1 is chosen for normalization, since this is 
the lightest one in our set of selected nuclei which has the two-parameter Fermi density 
distributions of nucleons. Being defined in this way, the transparency ratio better represents 
the systematic mass dependence for heavy nuclei. 

The transparency ratio Ra reveals strong local variations as a function of the mass number 
when calculated with the empirical nucleon density parameters (upper panel of Fig. H]). This 
arises from the details of empirical density profiles. For example, the local maximum for 
the 181 Ta nucleus appears due to the large diffuseness parameter of the charge distribution, 
a c h = 0.64 fm [20] . These local variations, as expected, disappear if we enforce the density 
profiles to be determined by the uniform parameters (lower panel of Fig. |4]). 

Another peculiar feature observed in Fig. @] (upper panel) is a strong drop of the trans- 
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FIG. 3. (Color online) The J/^f production cross section in p collisions with 27 A1, 5e Fe, 142 Ce 



56 t 



and Pb vs antiproton beam momentum calculated with ctj/^at = (dotted line) ctj/^jv = 3.5 
mb (solid line) and crj/^,^ = 6 mb (dashed line). Vertical arrows show the beam momentum 4.07 
GeV/c of the on-shell J production in vacuum. 



parency ratio along the isotope chain Sn. If we turn off absorption of both p and 

J/\l/, then the cross section for J/\l/ production at the on-shell peak varies along this isotope 
chain by about 5% only, because the proton density distribution is similar for the different 
isotopes. Therefore, this drop is mostly caused by p absorption and J/\l/ dissociation on the 
neutron excess in heavier isotopes. 

The sensitivity of the transparency ratio Ra to the input J / dissociation cross section 
is clearly visible in Fig. HI The only possible strong interaction channels of the J/\l/ dissoci- 
ation on a nucleon below DD production threshold (p t hr = 5.18 GeV/c) are J/^/N — >• A C D 
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FIG. 4. (Color online) The transparency ratio of J/^f production in antiproton-induced reactions 
gE) for the nuclei 9 Be, 12 C, 16 0, 27 Al, 40 Ca, 56 Fe, 63 Cu, 75 As, 112,116,120,124^ i42 Ce; i8l T&j i97 Au 
and 208 Pb plotted versus the mass number of the target nucleus. The ratio is normalized on 1 for 
27 A1. Open squares, solid and open circles represent calculations with (Jj/^n = 0, 3.5 mb, and 6 
mb, respectively. Upper panel - results with density parameters R q ,a q (q = n,p) determined for 
each nucleus separately as described at the beginning of sec. Lower panel - results for heavy 
nuclei excluding H2,ii6,i24g n with Rn = Rp = 1.16A 1 / 3 f m and a n = a p = 0.48 fm. 
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FIG. 5. (Color online) The ratio of the A C D production cross section to the J/vt production cross 
section at pi^b = 4.07 GeV/c for aj/^,^ = 3.5 mb (full blue circles) and for ctj/^tv — 6 mb (open 
purple circles). Heavy nuclei excluding 112 > 116 ,i24g n are s ] lown us i n g R n = R p = 1.16^4 1//3 fm and 
a n = a p = 0.48 fm. 

+ up to three pions. Hence, the J/^B dissociation cross section at the beam momentum of 
~ 4 GeV/c is equal to the inclusive A C D production cross section on the nucleon. More- 
over (see also 7|), both A c and D can not be absorbed in a nucleus for p\^ ~ 4 GeV/c 
but can only change momenta by rescattering on nucleons. For D, the absorption is ac- 
tually forbidden for any momentum by charm conservation in strong interactions. For A c , 
the threshold momentum in the nucleon rest frame for the A C N — > NND is 3.55 GeV/c, 
while the maximum A c momentum in J/tyN — > A C D is 3.34 GeV/c. The direct channels 
pp — > DD (p thr = 6.45 GeV/c) and pp —> A C A C (p thr = 10.16 GeV/c) are not reachable at 
the J/ty production threshold. Therefore, the cross section of the A c D-pair production in 
pA collisions at pi a b — 4 GeV/c can be simply calculated as 



o~a c d 



a 



w/oj/^abs 



pA->J/*(A-l)* a pA-^J/^(A-l)* > 



i w/o J/^abs. i ,i rn 

and (TpA-+j/q,(A-\y ~ b y the same Ec l- 



(26) 



), but with 



where a pA -+j/y(A-iy is given by Eq. 
-Pj/*,surv = 1- In Fig. [5] we show the ratio a Ac £,/apA^j/^(A-i)* at the on-shell peak of the 
J production vs target mass number. One sees the strong sensitivity of this ratio to the 
assumed value of the J/^/N dissociation cross section. 
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FIG. 6. (Color online) Transparency ratio Sa for the J/^f production in 7-induced reactions on 
nuclei vs target mass number. Left panel - = 20 GeV. Right panel - = 120 GeV. The results 
with crj/ifAT = 3.5 mb and oj/mN = 6 mb are shown by the solid and dashed lines, respectively. 
The lower and upper lines correspond to the calculations with formation length Ij/q, = and 
lj^ = 2 (12) fm for E 1 = 20 (120) GeV. Experimental data points from SLAC at £ 7 = 20 GeV 
(fBe and 181 Ta targets) Q and from Fermilab at £ 7 = 120 GeV (p, 9 Be, 56 Fe, and 208 Pb targets) 



111 ] represent the incoherent J/^> photoproduction cross section per nucleon normalized on 1 for 



9 



Be. 

We recall that the formation time effects are almost negligible and do not create an 
additional ambiguity for the J production in low-energetic antiproton-nucleus reactions. 
In contrast, formation time effects are very important for the 7-induced J/W production on 
nuclei. The transparency ratio in 7-induced reactions is denned according to Eq.flT]), which 
in the simplest approximation is expressed as (c.f. j^j): 

00 00 

Sa = <t^j,*x = 2tt 1 du I dzp ( Z:b)Pj/%smY ^ b) (27) 

^C r 7p->J/*X A- J J 
" ' —00 

This expression, however, is valid only at low photon energies, i.e. when the coherence length 
l c = 2£? 7 /m^ is much less than the nuclear radius. For E y = r 20 GeV and 120 GeV, where 



the J/ty production in photon- induced reactions is measured 10L 111], the coherence length 



is already quite large, l c = 0.8 fm and 4.8 fm, respectively. The deviations from the classical 
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probabilistic formula (12?]) appear in the Glauber model due to the coherent addition o 



f the 



production amplitudes on the two nucleons separated by the distance less than l c 
Refs. 



23 



23 



2J|. In 



241 ] , the formulas have been derived which generalize Eq. (12 7p for arbitrary values 



of l c (c.f. Eqs. (4.2a), (4.2b) in Ref.[23[ and Eq. (13) in Ref. (24j). The similar expressions 
are also given in Refs. 25j, |26j. Although these expressions somewhat differ from each other 
(mainly because of the different assumptions on the vector meson - nucleon elastic cross 
section), they all give the same limits of l c — > and l c — > oo. In the former case one gets 
Eq. (l27|) . while in the latter case one has to replace Pj/^ tSUTV (z, b) — > Pj/y )SU rv(— oo, b) in (|27[) . 
We will, thus, interpolate between these two limits by simply replacing Pj/y tSUTV (z, b) — > 
Pj/y,surv( z — h,b) in Eq. (l2TI) . This should be a quite rough approximation, but it serves at 
least our purposes of the exploratory studies of the photoproduction. 

Fig.[6]shows the mass dependence of the transparency ratio Sa corrected for the coherence 
length effects. We present results for the two previous values of the J/^N dissociation 
cross section, CTj/^at = 3.5 mb and 6 mb. However, the charmonium formation length 
= 2 (12) fm at pj/y = 20 (120) GeV/c is comparable with the nuclear size. Thus, the 
effective cross section 0"jyvj,^y of Eq. (1 1 1 1) is now substantially reduced with respect to cj/^j^ 
for the longitudinal coordinate within the nuclear target bulk region. The uncertainty in 
the determination of the charmonium formation length has an immediate feedback on the 
extraction of the genuine charmonium-nucleon dissociation cross section. As demonstrated in 
Fig. El for E 1 = 20 GeV, it is still possible to clearly see the difference between transparency 
ratios calculated with different values of ctj/^at. However, for E 1 = 120 GeV, the large 
formation length washes-out the sensitivity of Sa to ctj/vn- Moreover, the experimental 
errors do not allow to set tight constraints on Oj/qiN- 

Figures [7J and [S] present the beam momentum and mass number dependence of ^'(2S) 
production in p-induced reactions. Both dependences are quite similar to those for the J/^> 
production (c.f. Figs. |3]and|l]). The local variations of the transparency ratio (upper panel 
in Fig. [8]) due to the empirical density profiles are again visible. The smooth behavior of 
the transparency ratio as a function of the mass number is recovered if we substitute the 
empirical density parameters by the uniform ones (lower panel in Fig. [8]). 

To provide some hints on the possible charmonium absorption effects, we show in Figs. [7] 



and[8]the calculations with cr^'N = 20 mb as theoretically estimated in Ref. ^]. This reduces 
the yield by about a factor of 2-3 with respect to the calculation without absorption. 
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FIG. 7. (Color online) ^f'(2S) production cross section for p collisions with 27 A1, 56 Fe, 142 Ce and 
208 Pb vs antiproton beam momentum calculated with cr^/j^ = (dash-dotted line), c^/jv = 20 mb 
4^ with formation length l^i defined according to Eq. (jl3p (solid line) and with lq,> = (dotted 
line). Vertical arrows show the antiproton beam momentum of the on-shell ^'(2S) production in 
vacuum, p^b = 6.23 GeV. 

It is interesting to note, that with such a strong absorption, the production becomes 
almost independent on the target mass number. 

Since the beam momentum is now larger than for J /\& production, the formation length 
effects become visible (cf. solid and dotted lines in Fig. [7]). However, they are still weak as 
compared to the uncertainty caused by the largely unknown ty'N cross section. 

Finally, we address the multistep processes which are not included in the Glauber model. 



For this purpose we have performed the GiBUU model 12| calculations of the p+ 27 Al and 
p+ 208 Pb reactions. The GiBUU model takes into account the annihilation as well as the 
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FIG. 8. (Color online) The transparency ratio of VP' production in antiproton-induced reactions 
(I25I) as a function of the target mass number. Open squares show the calculation without 
absorption. Full and open circles represent the results with cr^'N = 20 mb [4| with ^'-formation 
length given by Eq. (fl~3j) at pi^ = 6.23 GeV/c and with l#i = 0, respectively. Upper panel 
- calculations with empirical density parameters R q ,a q (q = n,p) (see the beginning of sec. [3]). 
Lower panel - calculations for heavy nuclei excluding 112 > 11 M24g n uniform density parameters 
R n = R p = 1.16^ 1/3 fm and a n = a p = 0.48 fm. 
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FIG. 9. Comparison of the Glauber and GiBUU model calculations for the J/\I/ production cross 
section in p collisions with 27 A1 and 208 Pb vs beam momentum. Dashed lines correspond to our 
standard Glauber model calculations with the fixed proton energy E p = m — B as described in 
sect. [21 Solid lines are obtained within the Glauber model assuming the local Fermi gas for the 
nuclear ground state, same as in GiBUU. J/^f absorption cross section is set to zero. 



elastic and inelastic rescattering of the incoming antiproton with empirical cross sections. 
The nuclear density profiles are chosen to be identical in both, Glauber and GiBUU cal- 
culations. The antiproton-nucleon total and elastic cross sections fl91). f[T01) coincide in the 
both models as well. A comparison of the GiBUU and Glauber calculations is presented 
in Fig. [9j For simplicity, the nucleus was modeled in the local Fermi gas approximation in 
GiBUU. Therefore, for the comparison purposes we have also performed the Glauber model 
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calculations by doing the same assumption (solid lines in Fig. [9]). This has been achieved 
by replacing Eq. (]18p for the charmonium production width by the following formula: 



r 1 

1 p^R 



min(p 2 ,PF, P ) 2 + m 2 - JminQoi,^) 2 + m 



(28) 



Ap x?Lh E p q R 

where pi^ = \p\^(m 2 R — 2m 2 )^2EpmRqR\/2m 2 . The Glauber calculation is quite close to the 
GiBUU results at the peak. However, at higher beam momenta the Glauber model (with the 
local Fermi gas assumption) underpredicts GiBUU somewhat. The reason is that the fast 
antiproton has a chance to be decelerated by elastic or inelastic collisions with nucleons and 
get momentum closer to the peak momentum, where the cross section of J/^f production is 
larger (see also discussion in jg]). This mechanism is taken into account in GiBUU while it 
is neglected in the Glauber model. 



4. CONCLUSIONS 



We have performed the Glauber model calculations of J/fy(lS) and ^>'(2S) charmonium 
production in p-nucleus collisions at pi a b = 3 — 10 GeV/c. For the both charmonia, we have 
focused on the beam momentum range near the corresponding on-shell production peaks. 
Thus, only the pp — > R channel was taken into account. The main nuclear effect is the 
broadening and reduction of the narrow charmonium production peak due to the nuclear 
Fermi motion. 

The J/ty production cross section in pA collisions strongly depends on the input J/^/N 
dissociation cross section. This dependence is not blurred by the charmonium formation 
length, in contrast to J production in 7-induced reactions at high energies. 

The surface-dominated antiproton absorption leads to relatively large local variations of 
the J/^B transparency ratios as a function of the target mass number. This is due to the 
delicate interplay between the neutron and proton density profiles: The former governs the 
absorption range of the antiproton, while the latter defines the space region, where the J 
is produced. We conclude, that the quantitative determination of the J/^/N dissociation 
cross section from experimental data on J /^f production in pA reactions relies on the detailed 
and realistic description of the neutron and proton density distributions. 

We would like to recall at this point that the spreading of the proton momentum distri- 
bution due to the short-range correlations has not been taken into account in the present 
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calculations. Although these effects are very important at extreme kinematics regions, they 
will not sensitively modify our results near the on-shell peaks of charmonium production. 
The reason is that in this case the momentum integration in Eq.Q is not restricted from 
the low-momentum side, and, hence, is only weakly sensitive to the high-momentum tail. 
Another reason is that the short-range correlations become weaker at the nuclear surface, 
where the antiproton is predominantly absorbed. Overall, the short-range correlations may 
create the additional uncertainty of ~ 10% in the cross sections of the charmonium produc- 
tion close to the peak value. The same is true for the multistep effects due to the rescattering 
of the incoming antiproton in heavy target nuclei (Fig. 

It is expected from QCD, that the ty'N cross section is a factor of 2-4 larger than the 
J/^N cross section due to the larger size of the \P' as compared to the J and it may reach 



up to 20 mb 



271 ] . Such a strong absorption will be certainly testable with the new PANDA 



detector at FAIR starting from 2018. Having all the above uncertainties of our calculations 
in mind, we conclude that the measurements of the J/ty transparency ratio with a precision 
of at least ~ 20% would allow determination of the J/\l/-nucleon dissociation cross section 
with accuracy of about 3 mb. 

Significant J/^/N cross section implies the corresponding enhancement of the A C D pro- 
duction, since near J/W production threshold J/^/N — > A C D is the only possible inelastic 



channel of J/^/N interaction. On the other hand, there are severa 



spread predictions for the J/^/N — > A C D cross sections 22 
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models which give widely 



30j. Thus PANDA offers 



an interesting possibility to test these predictions by measuring the ratio of the A C D- to 
J/\l/-production cross sections (c.f. Fig. [5]). 
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